#include <iostream>
using namespace std;

double * AllocDouble1D(int n);